Upper-hybrid wave driven Alfvenic turbulence in magnetized dusty plasmas 
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The nonlinear dynamics of coupled electrostatic upper-hybrid (UH) and Alfven waves (AWs) is 
revisited in a magnetized electron-ion plasma with charged dust impurities. A pair of nonlinear 
equations that describe the interaction of UH wave envelopes (including the relativistic electron 
mass increase) and the density as well as the compressional magnetic field perturbations associated 
with the AWs is solved numerically to show that many coherent solitary patterns can be excited 
and saturated due to modulational instability of unstable UH waves. The evolution of these soli- 
tary patterns is also shown to appear in the states of spatiotemporal coherence, temporal as well 
as spatiotemporal chaos due to collision and fusion among the patterns in stochastic motion. Fur- 
thermore, these spatiotemporal features are demonstrated by the analysis of wavelet power spectra. 
It is found that a redistribution of wave energy takes place to higher harmonic modes with small 
wavelengths which, in turn, results into the onset of Alfvenic turbulence in dusty magnetoplasmas. 
Such a scenario can occur in the vicinity of Saturn's magnetosphere as many electrostatic solitary 
structures have been observed there by the Cassini spacecraft. 

PACS numbers: 52.35.-g, 52.25.Gj, 52.25.Vy, 52.25.Xz 



The nonlinear interaction of high-frequency (hf) elec- 
trostatic upper-hybrid (UH) waves (ES-UHWs) and the 
low-frequency (If) electrostatic or electromagnetic waves 
in a magnetized dusty plasma is a topic of current re- 
search (see, e.g., Refs. \^-§^) as those give rise to a great 
variety of nonlinear effects including parametric instabil- 
ities [7| , decay [S', 'S^ as well as modulational interactions 
[lji6^. The latter can be responsible for the evolution 
of UH envelope solitons [U |3l [5] and UH wave collapses 
[lOl[TT] which have been observed in laboratory and space 
plasmas [T0l[T2]. Moreover, the formation of such enve- 
lope solitons through the nonlinear interaction of coupled 
hf and If waves is one of the most interesting features in 
modern plasma physics in the context of plasma turbu- 
lence [IS'-TT, plasma heating and particle acceleration 
[18 . Furthermore, the pattern formation and the exis- 
tence of spatiotemporal chaos (STC) characterized by its 
extensive and irregular pattern dynamics in both space 
and time in nonlinear dynamical systems have received 
renewed interests (see, e.g. Refs. J3l HH |16j [23] ) . 

When the electric field intensity becomes strong and 
approaches the modified decay instability threshold, the 
interaction between hf and If waves results to 'weak tur- 
bulence' in which hf waves are scattered off If waves. 
However, if the electric field is so strong that it exceeds 
the modulational instability (MI) threshold, the interac- 
tion is then said to be in 'strong turbulence' regime in 
which transfer or redistribution of wave energy to higher 
harmonic modes with small wavelengths can take place 
[l3l[T4l[TC1. In this context, it is, however, believed that 
some sort of chaotic process may be responsible for the 
transfer of energy from large to small spatial scales [17j . 
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Recently, Shukla et al. [T have developed the non- 
linear theory of coupled mode interactions between ES- 
UHWs and modified Alfven waves (MAWs), and pre- 
dicted its application to the dusty magnetosphere of Sat- 
urn. However, they extended this work by accounting 
for the effects of the relativistic electron mass increase in 
the UH fields and the density as well as compressional 
magnetic field perturbations driven by the UH pondero- 
motive force. They made observations on a new class of 
oscillatory MI and spiky UH wave envelopes trapped in 
the electron density cavity. 

On the other hand, Williams et al [19] have reported 
observations of electric field solitary structures by Cassini 
spacecraft in the vicinity of Saturn's magnetosphere with 
ambient magnetic fields ~ 0.1 to 8000 nT. They also re- 
ported a range of the peak-to-peak wave amplitude by the 
external magnetic field as well as time durations of the 
solitary pulses. Moreover, it has been suggested that the 
turbulence in the Earth's magnetosheath may account 
for the observations of many solitary pulses there [20] . 
Also, majority of solitary structures have been found at 
Europa in the moon's wake, a region of known turbulence 
[2T] . It is, therefore, of natural interest to investigate the 
possibility for the occurrence of turbulence in the vicinity 
of dusty magnetosphere of Saturn where many ES-UHWs 
have been observed [19 . We report here, possibly for the 
first time, the onset of If Alfveenic turbulence in the non- 
linear interaction of ES-UHWs and MAWs in parameter 
regimes that are relevant to dusty magnetosphere of Sat- 
urn. We find that redistribution of wave energy indeed 
takes place to higher harmonic modes, and the process 
becomes faster the smaller is the wave number of mod- 
ulation at which many unstable modes are to be excited 
and saturated by the MI of ES-UHWs. 

In what follows, we consider the following nondimen- 
sional set of equations describing the dynamics of coupled 
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ES-UHWs and MAWs [J. 

i [dt + aid^) E + a2dlE + ag ( 
{d^,+h,-hdl)N-hdl\E\^ 
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N)E = 0, 
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where the coefficients are ai = ZhQaVre/^H^ cl2 = 

= m/ic^/3F^g in which ko is 
= 1/ (l - Scj^/cj^) > with ujc{p) de- 
noting the electron cyclotron (plasma) frequency, Vtc is 

the electron thermal speed, ujh = 'sj^p + UH 

resonance frequency and uOpi is the ion plasma frequency. 

Also, VLq = 'sjojp + 2uj1 , r^i^ = ZdfidoUJci/'^eO IS the Rao 

cut-off frequency [22 , is the number of electrons re- 
siding on a stationary charged dust grain with unper- 
turbed number density Udo, Var = /i^o/ V^Tm^o^ is 
the modified Alfven speed [22], /i = nio/^eo(= I -\- S = 
1 + ZdTido/neo) is the unperturbed ion to electron num- 
ber density ratio, m = rrie/mi is the electron rest to ion 
mass ratio and Bq is the external magnetic field. More- 
over, E = Ex/Ec is the x-component of the nondimen- 
sional UH wave electric field with Ec = 2meCOJp/^/3eujH 
and TV = nei/neo(= Biz/Bq <C 1) denoting the nor- 
malized electron number density perturbation (compres- 
sional magnetic field) associated with the MAWs in dusty 
plasmas [5 . The space and time variables are normalized 
respectively by Vre/^pi and the ion plasma period a;"/. 
The fourth and fifth term oc in Eq. arise respec- 
tively due to the relativistic electron mass increase and 
the density as well as the compressional magnetic field 
fluctuations that are driven by the ES-UHW pondero- 
motive force [For details about Eqs. ([T]) and ^ see, e.g. 
Ref. [1]]. 

Now, the dispersion relation for the MI of 
a constant amplitude UH pump of the form 
^0 exp {ikox — ico'i:/t)+complex conjugate, can readily be 
obtained. Here we decompose the electric field E as the 
sum of the pump and its two sidebands, and obtain the 
following dispersion law pLj (Note that some terms in Ref. 
[1 might be missing or have been simplified somehow). 
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2a^El),K = a2K^ 



= 2asbsK'^E^A, where 
-asEln\^ = h^b2K^^ 



Vg = ai are in nondimensional forms and 1], K are the 
normalized wave frequency and wave number of modu- 
lation respectively. For the parameter values as relevant 
to the Saturn's magnet osphere, namely ngo = 80 cm~^, 
^ = 10, Te = 3 X 10^ K, Bo = 1500 nT [19] we find that 
the coefficients of and Q become smaller compared 
to the other terms. Hence the above dispersion relation 
can be simplied to obtain ^'^ ^ ^ [^ar - K'^Vg + A] ± 

\ ( + K^V^ - A] V ^a,h,K^ElK] 



This 



shows that the ES-UHWs are modulationally unstable 
when the expression inside the square root beco mes ne g- 
ative (except the case in which A > for > ^/a^Ja2E{) 
giving a purely damping mode for the lower sign) for 
K < Kc^ a critical wave number. We numerically 




FIG. 1. (Color online) Contour of \E[x, t)\ =const for K = 1.1 
showing that the pattern selection leads to only one harmonic 
pattern. The syetm is in the coexistence of SPC and TC. 
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FIG. 2. (Color online) Contour of \E{x,t)\ =const for K = 
0.52. Four solitary patterns collide pairwise at t ~ 45 and 
t 90 to fuse into two new incoherent patterns. The syetm 
is still in the coexistence of TC and STC. 



investigate some domains of K with the parameters as 
above for different values of the initial pump for which 
the MI sets in. These are namely, 0.23 ^ K ^ 0.51 for 
^0 = 0.2; 0.47 <K< 0.83 for = 0.3; 0.74 <K<l.l 
for £^0 = 0.4; 1.1 < < 1.4 for Eq = 0.5 etc. In these 
regimes only one solitary pattern may be formed by the 
master mode due to MI, and beyond which many modes 
will be excited and saturated by the wave number as 
can be seen from following simulation results. 

We numerically solve the Eqs. ([T]) and ([2| by a fourth- 
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FIG. 3. (Color online) Contour of \E(x,t)\ =const for K = 
0.1. Eighteen solitary patterns are initially formed from the 
pattern slection. After sometime, they collide and fuse to 
form only six new incoherent patterns. The Alfven wave emis- 
sion also occurs everywhere. The collision is random and not 
confined between two patterns. The STC state emerges. 



order Runge-Kutta scheme where the spatial derivatives 
are approximated with centered second-order diference 
approximations. In the numerical scheme we consider 
the time step dt = 10~^ and assume spatial periodicity 
with the simulation box length = 27t/K (the res- 
onant wavelength) together with the grid size 1024 or 
2048 depending on so that x = corresponds to the 
grid position 512 or 1024. We choose the initial condi- 
tions as [El [H [16] E{x,0) = Eo^ Ei cos (Kx) , N = 
Ni cos {Kx), where Ei = b{bi^b2K'^), Ni = -2bb3EoK'^ 
and b = 1/500 is an arbitrary constant to ensure that the 
perturbation is small. Basically, for Kc/2 < K < Kc-, one 
unstable mode with spatially modulational scale = 
27t/K is excited and then saturated to form few spatially 
periodic solitary patterns. However, as K is lowered from 
there may exist many solitary patterns with spatial 
length scales Im = Lx/m, where m = 1 is for the mas- 
ter mode and m = 2, 3, M for the unstable harmonic 
modes with M < Mp = [Kc/K] being due to the pat- 
tern selection. Thus, the envelope E can be expressed 
as [H m Hi [23] ^ = '£^^^Em{t)exp{imKx) + 
Em{t) ex.-p {imKx), in which the first summa- 
tion comes from the master mode and the unstable har- 
monic modes, whereas the second one is due to the non- 
linear interaction of hf and If modes. We now proceed 
with the numerical simulation of Eqs. ([T]) and ^ con- 
sidering Eq = 0.5 for which the domain of MI for K is 
1-1 ^ ^ ^ 1-4. In the region Kc/2 = 0.7 < K < 1.4, the 
motion of the solitary pattern is either temporal recurrent 
(periodic) or pseudorecurrent (quasiperiodic) depending 



on the values of K. 

Figure 1 shows that for K = 1.1, only one master 
pattern is formed ai x ^ 180 implying that the system 
is stable. However, as can be seen from Fig. 2 that 
if K is lowered, i.e., for K = 0.52 < 0.7, one master 
pattern and one harmonic pattern whose initial peaks are 
dit X ^ 750 and 550 exhibit stochastic behaviors, which 
results in the coexistence of temporal chaos (TC) and 
spatially partial coherence (SPG). It has been shown that 
the resonant overlapping may be the cause for TG [23] , 
Furthermore, in Fig. 2, the pattern selection leads to the 
excitation of four solitary modes formed initially from 
the master mode and unstable harmonic modes. At t ^ 
45, two solitary patterns initially peaked at x 50 and 
350 collide and fuse to form a new pattern with stronger 
amplitude. At the same time strong Alfven wave (AW) 
emission takes place everywhere and solitary structures 
are distorted. As a result, at t 90, one solitary pattern 
peaked ai x ^ 550 collides with the master pattern and 
gets fused into another new one. Thus, after the two 
collisions there remains only one new incoherent pattern 
together with the distorted master pattern. However, as 
will be evident from the analysis of power spectra (Figs. 
4 to 6 below) that few collisions are not sufficient for 
the cause of STG. The coherence of the system is still 
partially retained, so that the system is in the state of 
TG and SPG. 

Next, we consider the case in which many unstable 
modes are excited and saturated from the master mode, 
unstable harmonic modes as well as due to nonlinear in- 
teractions to form many solitary patterns. As for exam- 
ple, for K = 0.1, the pattern slection leads to eighteen 
solitary modes (Fig. 3) in which the fisrt collision occurs 
at t ~ 45. At longer times only few new incoherent pat- 
terns remain after collision and fusion among them due 
to strong AW emission. Four patterns initially peaked at 
X ^ 1700, 1800, 1900, 2000 collide pairwise dit t ^ 6b 
and 80 respectively to form two new patterns, which 
again collide to each other due to strong AW emission 
at t 130 to fuse into another new pattern. Also, the 
harmonic pattern which was initially excited at x ~ 1200 
disappears after some time t ^ 110 due to AW emission. 
As time goes on, several other collisions and fusions take 
place repeatedly. The patterns are then much distorted 
and the original eighteen solitary waves are finally fused 
into six new incoherent patterns. 

We now perform a wavelet analysis of the electric field 
time series using a Morlet transform. Wavelets have been 
widely used for the analysis of time-frequency spectra 
in diverse applications including dynamical analysis of 
structural systems, image processing as well as pattern 
recognition in spatiotemporal systems [24 . In order to 
recognize the spatiotemporal features of the present sys- 
tem as observed in Figs. 1 to 3, we analyze the corre- 
sponding data of the electric field for different values of 
K = 1.1, 0.52 and 0.1 as above. To this end, we rescale 
the absolute value of the electric field amplitude as {Ew) 
by the standard deviation of the sampling data and plot 
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FIG. 4. (Color online) The normalized wave electric field (Ew, 
upper panel) rescaled with the standard deviation of the sam- 
pling data and the wavelet power spectra {Pw in logarithmic 
scale, lower panel) with respect to the sampling time T for 
K = 1.1 indicating periodic wave trains. 



against the sampling time T (with a time step 0.01) which 
depends on the size of the data considered. In our analy- 
sis we have used a Morlet wavelet as the mother wavelet 
which offers a good balance between time and frequency 
localizations and consists of a plane wave modulated by 
a Gaussian wave function. The Morlet transform is then 
used to calculate the wavelet for the field at the signif- 
icance level with variance, = 1. The corresponding 
wavelet power spectra of the time series shown in the 
upper panels of Figs. 4 to 6 are visualized in the lower 
panels of the same. We see that in a certain interval of 
T, the pattern shows a regular simple periodic train of 
waves for if = 1.1 (Fig. 4). This result is consistent with 
the contour plot as shown in Fig. 1. When K = 0.52, 
Fig. 5 indicates a quasiperiodic wave trains, i.e., chaotic 
in time but may be partially periodic in space as analo- 
gous to the behaviors depicted in Fig. 2. Finally, strong 
chaotic motion can be seen from Fig. 6 for = 0.1 which 
agrees with the features observed in Fig. 3. 

Thus, in the process of collision and fusion strong 
chaotic activity in both space and time can be observed 
in the present system. As a result, there remains, af- 
ter a longer time, few incoherent patterns together with 
few stable higher harmonic modes excited through non- 
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FIG. 5. (Color online) The same as in Fig. 4, but for K = 0.52 
indicating the state of temporal chaos but spatially partial 
coherence. 



linear interactions. A certain amount of the system en- 
ergy, which was initially confined to many solitary waves, 
is now redistributed into those new incoherent patterns 
as well as to stable higher-harmonic modes with short- 
wavelengths {K > Kc). The transport of the trapped 
electrons then occurs, and electrons can thus be acceler- 
ated across the magnetic field via a surfatron acceleration 
mechanism [25]. So, if initially there exist a number of 
unstable modulation lengths to form many solitary pat- 
terns, pattern evolution (i.e., collision, fusion and distor- 
tion among them) and strong AW emission can lead to 
the STC state and hence to the onset of turbulence in 
dusty magnetoplasmas. However, the mechanism that 
leads to the state of STC is still remains unclear. 

We ought to mention that in the present investiga- 
tion, we have considered the parameter regimes that give 
positive group dispersion {a > 0) of the UH waves. In 
the opposite case, i.e. for a < 0, one needs to consider 
some higher magnetic field strength than that considered 
here. This may, however, give rise stable wave propaga- 
tion |3[26|. Furthermore, the effect of the cubic nonlin- 
ear term in Eq. ([T]), which appears due to the relativis- 
tic electron mass increase, is to enhance the excitation 
of many more unstable harmonic modes and hence the 
faster dynamical transition from order to STC, than the 
nonrelativistic case On the other hand, the stationary 
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FIG. 6. (Color online) The same as in Fig. 4, but for K = 0.1 
indicating that the motion is strongly chaotic in both space 
and time. 



charged dust grain, which modifies the frequency of the 



modulated ES-UHWs and hence the associated instabil- 
ity growth rate, does not change the qualitative behaviors 
(except the fact that the collision and fusion among the 
patterns somewhat differ at different times and that the 
wave amplitude decreases slightly at a specific time) of 
the nonlinear interactions as shown in Fig. 3. That is, 
the number of modes initially excited at different modu- 
lational length scales and the number of remaining modes 
(i.e., six) after collision and fusion among them, remain 
the same. Moreover, the nonlinear coupling of collid- 
ing multiple UH wave envelopes (head-on collisions) that 
may give rise to stability of one-dimensional envelope 
solitons due to negative group dispersion (a < 0), could 
be an another important investigation [26] , but is beyond 
the scope of the present study. 

It is to be mentioned that although, there is no di- 
rect or indirect observations, in particular, for localized 
UH wave envelopes coupled with Alfven waves. However, 
based on the recent observations, i.e. existence of many 
electric field solitary structures in the vicinity of Saturn 
magnetosphere with the magnetic fields 0.1 nT to 8000 
nT [19 , and also to the fact that the equatorial plane 
of Saturn contains charged dusts, and dusty plasma can 
move across the magnetic field lines [2 , we conclude that 
the present results could give an insight for understand- 
ing the origin of If Alfvenic turbulence that can occur in 
Saturn's magnetosphere. However, as of 2010, the mag- 
netosphere of Saturn still remains an important subject 
of the ongoing investigation by Cassini mission. 
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